TU-812 
April, 2008 



Big-Bang Nucleosynthesis and Gravitino 

o 

(N 

u 

Oh 

( a )( 6 )Masahiro Kawasaki, (^Kazunori Kohri, ^Takeo Moroi, 

r^\ ' ill 

C\| ! ^Akira Yotsuyanagi 

^Institute for Cosmic Ray Research, University of Tokyo, Kashiwa 277-8582, JAPAN 

^Institute for the Physics and Mathematics of the Universe, University of Tokyo, 

Kashiwa 277-8582, JAPAN 

j> ■ ^Physics Department, Lancaster University, Lancaster LAI J^YB, UK 

in 

^ \ ^Department of Physics, Tohoku University, Sendai 980-8578, JAPAN 

m 

o 

oo 
o 



X 



Abstract 

We derive big-bang nucleosynthesis (BBN) constraints on both unstable and stable 
gravitino taking account of recent progresses in theoretical study of the BBN processes 
as well as observations of primordial light-element abundances. In the case of un- 
stable gravitino, we set the upper limit on the reheating temperature assuming that 
the primordial gravitinos are mainly produced by the scattering processes of thermal 
particles. For stable gravitino, we consider Bino, stau and sneutrino as the next-to-the- 
lightest supersymmetric particle and obtain constraints on their properties. Compared 
with the previous works, we improved the following points: (i) we use the most recent 
observational data, (ii) for gravitino production, we include contribution of the longi- 
tudinal component, and (iii) for the case with unstable long-lived stau, we estimate 
the bound-state effect of stau accurately by solving the Boltzmann equation. 



1 Introduction 

It is well known that, in the framework of the standard model of particle physics, viable 
scenario of the evolution of the universe cannot be obtained. One of the most important 
reasons is that, in the standard model, there is no candidate of dark matter. Among various 
possibilities, the lightest superparticle (LSP) in supersymmetric models is a good candi- 
date of dark matter, and hence supersymmetry is thought to be very attractive from the 
cosmological point of view as well as from particle-physics point of view. 

In considering cosmology based on supersymmetric models with i?-parity conservation, 
it is important to be aware that graviton also has its superpartner, which is called gravitino. 
Even though gravitino is an extremely weakly interacting particle, it may play a significant 
role in the evolution of the universe. In particular, with the existence of gravitino, there may 
exist a long-lived superparticle whose lifetime is longer than ~ 0.1 sec. For example, if one 
of the superparticles of standard-model particles is the LSP, gravitino becomes unstable and 
is very long-lived. On the contrary, if the gravitino is the LSP, lightest superparticle in the 
minimal supersymmetric standard model (MSSM) sector, which we call MSSM-LSP, decays 
into gravitino (and something else) with a Planck-scale suppressed decay width. 

Long-lived superparticles may affect the big-bang nucleosynthesis (BBN) if their life- 
times are longer than ~ 0.1 sec. Such a long-lived superparticle is produced in the early 
universe and decays in the thermal bath at the cosmic time comparable to its lifetime. 
Importantly, the energetic particles emitted by the decay may induce hadronic and electro- 
magnetic shower and produce large amount of energetic hadrons and/or photons. Those 
energetic hadrons and photons induce hadro- and photo-dissociation processes of light ele- 
ments which are synthesized by the standard BBN (SBBN) reactions. In addition, mesons 
such as charged pions are also produced in the hadronic shower induced by the decay of 
long-lived particle. Such charged pions and baryons cause extraordinary p <-> n conversion 
processes, resulting in the change of 4 He abundance. Since the predictions of the SBBN 
scenario are in reasonable agreements with observations, we obtain upper bound on the 
primordial abundance of the long-lived superparticle. 

In this paper, we study the BBN scenario in detail for the cases where gravitino is 
unstable and stable, and derive constraints on those cases. In supersymmetric model, such 
BBN constraints have been extensively studied for the cases where gravitino is unstable 
PISElHlElEjmE] and stable [101 [LTJ [[2j [13] . The primary purpose of this paper is 
to study the constraints on scenarios with long-lived superparticles from the BBN scenario, 
taking account of recent developments in this field. Compared to the most recent studies 
[SI El El EH EH] , we have improved the following points: 

• We have used the most recent observational constraints on the primordial abundances 
of light elements. 

• We have taken into account the contribution of the production of the longitudinal 
component of gravitino in the calculation of the primordial abundance of gravitino. 



Consequently, we have obtained severer upper bound on the reheating temperature for 
the case of unstable gravitino, in particular when the gravitino mass is relatively light. 

• For the case that the gravitino is the LSP and a charged slepton is the next-to-lightest 
super symmetric particle (NLSP), we have calculated the abundance of the bound state 
of the charged slepton with 4 He solving the Boltzmann equation without the use of 
Saha formula and hence have estimated the catalyzed production of 6 Li precisely. In 
addition, for the scenario with Bino-like NLSP, we have performed a detailed analysis. 

The organization of this paper is as follows. In Section [2, we first summarize how the 
BBN processes are analyzed taking into account the effects of long-lived superparticles. 
Then, in Section El we discuss constraints on the case where gravitino is unstable. We pay 
particular attention to the scenario where the LSP is the lightest neutralino and derive the 
upper bound on the reheating temperature after inflation. In Section [H we consider the case 
where gravitino is the LSP and hence is stable. In such a case, the MSSM-LSP becomes 
unstable and its decay may affect the BBN scenario. The constraints depends on what the 
MSSM-LSP is, so we consider several cases where the MSSM-LSP is the lightest neutralino, 
one of the charged slepton, and one of the sneutrino. Section [5] is devoted to conclusions and 
discussion. 

2 BBN with Long-Lived Unstable Particles 

2.1 Effects of the long-lived particle and calculation procedure 

We first summarize the general procedure to study effects of long-lived unstable particles on 
BBN. In particular, in this section, we explain how we calculate non-standard BBN reaction 
rates with long-lived unstable particles. We also summarize observational constraints on 
light-element abundances which will be adopted in our study. 

In our study, we consider two important types of decay modes of long-lived particles: ra- 
diative decay and hadronic decay. With the radiative decay, photon and/or charged particles 
are emitted. Then the light elements would be photo-disintegrated (or photo-dissociated) 
by energetic photons which are produced in the electro-magnetic shower induced by high en- 
ergy photons/charged particles [21 El HI El US] • With the hadronic decay, the emitted hadrons 
induce extraordinary interconversions between background proton and neutron (p <-► n con- 
version) through the strong interaction. Such conversion processes enhance the neutron-to- 
proton ratio, resulting in the over-production of 4 He [HI El El El El- The energetic nucleons 
also destroy the background He and nonthermally produce other light elements, such as D, 
T, 3 He, 6 Li, 7 Li and 7 Be [HI [191 El El El E] • Furthermore, if the long-lived particle is charged, 
it may form a bound state with background nuclei and change the nuclear reaction rates. 
In particular, for the production of 6 Li, such a process, which is called catalyzed production 
process [20], may become very important. 

We incorporate these effects, photo-dissociation, p <-> n conversion, hadro- dissociation, 
and catalyzed production process for the case with charged long-lived particle, into a stan- 



dard BBN code [21] and calculate light-element abundances. (Some details of these effects 
are discussed in the following.) In order to estimate errors in theoretical calculation, we 
perform Monte-Carlo simulation by changing both the standard and non-standard reaction 
rates within their statistical and systematic errors. (For the baryon-to-photon ratio, we use 
77 = 6.10 ± 0.21 j22j.) The estimated errors in the calculation of light-element abundances 
are used in the estimation of the \ 2 variable to derive BBN constraints; in the calculation 
of the x 2 variable, theoretical and observational errors are added in quadrature. 

2.1.1 Radiative decay 

In order to study the effects of radiative decay, we numerically solve the Boltzmann equation 
governing the time evolution of the photon and electron spectra. To obtain photo-dissociation 
rates, we take convolutions of the photon spectrum and the energy-dependent cross sections. 
In the following, we briefly discuss our treatment of the photo-dissociation processes. (For 
the details of the effects of radiative decay, see [3].) 

If the parent particle radiatively decays into high-energy photons or charged particles, 
those daughter particles scatter off the background photons 7bg and electrons e# G via 
electromagnetic interaction. (Here, the subscript "BG" is for background particles in the 
thermal bath.) The most important processes are the Compton scattering (7 + eg G — > 
7 + e~), the electron-positron pair creation (7 + 7bg — ► e~ + e + ), the pair creation in nuclei 
(7 + iV — » e~ + e + + N, with N being the background nuclei), the photon-photon scattering 
(7 + 7bg - > 7 + 7), an d the inverse Compton scattering (e~ + 7bg — > e~ +7). Consequently, 
the electromagnetic cascade shower is induced and a series of these processes produce many 
photons. 

The timescale of the energy loss processes of the photons is much faster than the timescale 
of cosmic expansion if the cosmic temperature is sufficiently high T > 0(1 eV). On the 
other hand, the thermalization timescale of the photons can be longer than that to react 
with the background light elements, and hence the photons may cause the photo-dissociation 
processes; the photo-dissociation process with the threshold energy of Q becomes effective 
when the cosmic temperature is lower than ~ m? e /22Q [4J. For the photo-dissociation of 
D ( 4 He), for which the threshold energy is Q ~ 2.2keV (20 MeV), the cosmic temperature 
should be lower than ~ lOkeV (~ IkeV). 

In addition, once the 4 He is destructed, energetic T and 3 He are produced. They scat- 
ter off background 4 He and produce 6 Li nonthermally. We also include photo-dissociation 
processes of 6 Li, 7 Li and 7 Be. 

2.1.2 Hadronic decay 

Next, we summarize how the effects of hadronic decay are studied. If the parent particle 
decays into quarks and/or gluons, they immediately fragment into high-energy hadrons. 
Comparing the lifetime of the produced hadrons with the timescale of their interactions 
with the background nuclei, only long-lived mesons, 7r =t , K^, K^, as well as the nucleons 



n and p are important for the study of the BBN reactions; hadrons with shorter lifetimes 
decay before scattering off the background nuclei. 

We calculate the distributions of partons emitted from the decaying parent particle. 
In particular, for the systematic study of the decay chain, we use the ISAJET/ISASUSY 
packages [23] . Then fragmentation processes of the partons into hadrons are simulated by 
using the event generator PYTHIA [21] and finally the distributions of the produced hadrons 
are obtained. 

The long-lived (and stable) hadrons destroy the background 4 He and produce D, T, 
and 3 He copiously. Furthermore, the energetic daughter T, 3 He and the scattered 4 He also 
nonthermally produce 6 Li, 7 Li and 7 Be through the collisions with the background 4 He. 
These processes are severely constrained by the observations [IH El E] . 

On the other hand, the neutron injection at around T ~ 30 keV may reduce the final 
abundance of 7 Be (or 7 Li) through the neutron capture process, 7 Be(n,p) 7 Li(p,a) 4 He [25| |6| 
Q2J ES]- However, this process does not affect our constraints because we adopt a very 
mild observational abundance of 7 Li. 

With the hadronic decay, the conversion processes of the background p and n also occur. 
In particular, once protons, neutrons, and pions (and their anti-particles) are produced in 
the hadronic shower, they induce p <-» n conversion processes and change the neutron to 
proton ratio n/p even after the normal freeze-out epoch of neutronej Such p «-► n conver- 
sion processes increase 4 He and hence are constrained from observation of the primordial 
abundance of 4 He. 

We solve the time-evolution of the hadron distribution functions including all the relevant 
electromagnetic and hadronic energy-loss processes, and calculate the reaction rates of the 
light elements taking into account the hadro-dissociation and p <-> n conversion processes 
(as well as those of the photo-dissociation). For further details, readers can refer to [B|. 

Before closing this subsection, we discuss the minor modification of the treatment of the 
hadronic shower compared to the previous studies [5j El [7J . In the present study, we adopted 
a new evolution scheme of the hadronic shower. For the study of the abundance of non- 
thermally produced 6 Li, 7 Li and 7 Be, it is important to obtain information on the energy 
distribution of 4 He produced by the scattering process N + 4 He — > N' + 4 He + 7r. However, 
there is a lack of experimental data of the energy distribution of the scattered 4 He in the 
inelastic N + 4 He collision. In the old papers [5J El [TJ, quantum chromodynamics (QCD) 
properties was used to extrapolate results of high-energy experiments to lower energies (as 
far as the extrapolation is kinematically allowed). Consequently, the final-state 4 He acquired 
larger energy than that estimated via equipartition distribution in the center-of-mass system. 
However, in the present study, we adopted the smaller value of the 4 He energy among 
the equipartition value and the QCD prediction. This procedure gives more conservative 
constraints. For non-relativistic iV + 4 He collisions, this modified method gives smaller 
energies of the scattered 4 He, which reduces the resultant nonthermally produced Li and Be 
by ~ 50 % or so [26]. However, this does not affect the constraints given in [5j El [7J. That 



^We omit contributions from K^ and K\ to this p <-> n conversion according to a discussion in 6 
although there have been some attempts to include effects of Kaons [121 [5] . 



is because the nonthermally produced Li and Be are subdominant compared to the SBBN 
contributions, and also because we adopt conservative observational constraints. 

2.2 Observational light-element abundances 

Next, we show observational constants on primordial abundances of D, 3 He 4 He, 6 Li and 7 Li, 
which we adopt in our study. The errors are presented at 1 a. The subscript "p" and "obs" 
are for the primordial and observational values, respectively. 

We adopt the following observational constraints on the deuterium abundance 

(n D /n H ) p = (2.82 ± 0.26) x 1(T 5 , (2.1) 

which is the most-recently reported value of the weighted mean [27J. This value well agrees 
with the baryon to photon ratio suggested by the WMAP 3-year CMB anisotropy observa- 
tion [22]. 

To constrain the primordial 3 He abundance, we use an observational 3 He to D ratio as 
an upper bound, which is a monotonically increasing function of the cosmic time. (For 
details, see [2HI El IB]-) In this study we adopt the newest values of D and 3 He abundances 
simultaneously observed in protosolar clouds (PSCs), (n3 He / , n.H)psc — (1-66 ± 0.06) x 10~ 5 
and (n D /n H )psc = (2.00 ± 0.35) x 10~ 5 [29J. Then we get 

{m He /n D ) p < 0.83 + 0.27. (2.2) 

Concerning the mass fraction of 4 He, in [30], two values are recently reported by using old 
and new 4 He-emissivity data, Y p = 0.2472 ± 0.0012 and Y p = 0.2516 ± 0.0011, respectively. 
It should be noted that the latter value is inconsistent with the SBBN prediction (~ 0.2484) 
even adopting possible theoretical errors of 0.0004. However, we do not take its face value 
since the error presented in [30] does not include systematic effects. Thus, in this study, we 
add an error of 0.0040 [31J to derive conservative constraint: 

F p = 0.2516 ± 0.0040. (2.3) 

(For the systematic uncertainty of the observed 4 He abundance, see also |32j.) 



For 7 Li, we adopt the most recent value of the 7 Li to hydrogen ratio log 10 ( 7 Li/H) 



obs 



—9.90 ± 0.09 given in [33J. This is close to the value given in [34]: Logi ( 7 Li/H) obs = 
—9.91 ± 0.10, and in addition, slightly larger values have been also reported for years; for 
example, see [35]: Logio( 7 Li/H) OD s = —9.63 ± 0.06. Those observational values are smaller 
than the SBBN prediction by approximately 0.3 dex or so. Concerning the inconsistency 
between the face value of the primordial 7 Li abundance and the SBBN prediction, various 
solutions have been discussed from the viewpoint of both astrophysics [36], [37] and cosmology 
[381 I2H1[71[121 H71 [391 HU SQl HH [261 SH SSlIlll. However, the main purpose of this paper is 
to derive a conservative constraint, so we do not go into the details of these models. Instead, 
we add an additional systematic error of +0.3 dex into the observational face value m Li /nn 
in our study: 

Log 10 (n7 Li /n H ) P = -9.90 ± 0.09 + 0.3. (2.4) 



Notice that +0.3 dex is the systematic error expected from effect of 7 Li depletion by rotational 



#2 



mixing in stars [36 

As for a 6 Li constraint, we use {nsix/mii^oha = 0.046 ± 0.022, which was newly-observed 
in a very metal-poor star [47] . We also add a systematic error of +0.084 [13] to take into 
account depletion effects in stars adopting the relation between 7 Li and 6 Li depletion factors, 
A log 10 (n6 Li /nH) = 2.5Alog 10 (n7Li/n H ) [H1I36], which leads to A log 10 (n6 Li /n7 U ) = 0.45 with 
7 Li depletion factor 0.3 dex. Then we get the following upper bound: 

{ne u /nr Li ) p < 0.046 ± 0.022 + 0.084. (2.5) 

3 Unstable Gravitino 

3.1 Primordial abundance of gravitino 

When gravitino is unstable, gravitinos produced in the early universe may affect the light- 
element abundances because the lifetime of gravitino is extremely long. Indeed, if gravitino is 
lighter than ~ O(10 TeV), its lifetime becomes longer than ~ lsec and, in such a case, most 
of the primordial gravitinos decay after the BBN starts. Since the predictions of the SBBN 
are in reasonable agreements with observations, the abundance of gravitino is required to be 
small enough. Thus, we obtain upper bound on the abundance of primordial gravitino. 

There are several possible mechanisms of producing gravitino in the early universe; grav- 
itino may be from the decay of condensations of scalar fields, like moduli fields [49] or inflaton 
field [50], as well as from scattering processes of particles in the thermal bath. Gravitino 
production due to the decay of scalar condensations is, however, highly model-dependent. In 
order to derive conservative constraints, we presume that the primordial gravitinos originate 
from the scattering processes. In such a case, the abundance of gravitino is approximately 
proportional to the reheating temperature after inflation. Then, in order not to spoil the 
success of BBN scenario, upper bound on the reheating temperature is derived. 

Concentrating on the gravitino production by the scattering, the evolution of the number 
density of gravitino is governed by the following Boltzmann equation: 

^ = -3Hn 3/2 + C 3/2 , (3.1) 

where H is the expansion rate of the universe and C3/2 is the collision term. In our analysis, 



we study the above Boltzmann equation simultaneously with other relevant equations^ 3 

dpinf 



dt 

dprad 

dt 



3-f/pinf — TinfPinf, (3.2) 

AHp raxi + r inf p inf , (3.3) 



# 2 See also [151135] for detailed discussion on possible systematic uncertainties. 

# 3 Here, we assume that the potential of the inflaton field is well approximated by parabolic potential. 

6 



where pi n f and p ra d are energy densities of inflation and radiation, respectively, and r in f is the 
decay rate of the inflaton field. We numerically solve Eqs. (13. ip — (13. 3p with initial conditions 
n z/2 = and p in f ~^> px', we follow the evolution of 713/2, Pint, and p ra d from the inflaton- 
dominated epoch (i.e., t <C r^ f ) to the radiation-dominated epoch (i.e., t 3> r^). Most of 
the energy density stored in the inflaton is converted to that of radiation at the reheating 
temperature which is of the order of y/M^T^ (with M* ~ 2.4 x 10 18 GeV being the reduced 
Planck scale); in this paper, the reheating temperature is defined as 3H(T R ) = T in{ which 
gives 

( 10 \ 1/4 

Tr= (—, 7f^- 2 MtTlA , (3.4) 

\9*{Tr)^ 2 J 

where ^*(Tr) is the effective number of the massless degrees of freedom at the time of 
reheating. In our study, we use <?*(T R ) = 228.75. 

The collision term has been calculated in [SH E2] by taking account of the thermal masses 
for gauge-boson propagators: 

^ISKSj?"*^ 1+ *47J ln ^(r)' (3 ' 5) 

where the summation is over the standard-model gauge groups, and c a = 11, 27, and 72, 
while k a = 1.266, 1.312, and 1.271, for U{1) Y , SU(2) L , and SU(3) C , respectively^ In 
addition, M a (a = 1 — 3) denotes the gaugino mass. Compared to the previous works on 
unstable gravitino [HI [7], we have included the production of the longitudinal mode, (part 
of) which is from the term proportional to M^/m\i 2 in Eq. ( 13. 5ft . Although we numerically 
calculate the primordial abundance of gravitino, we also provide its fitting formula. Let us 
define the yield variable for particle X with lifetime tx as 



Y x = 



n x 



(3.6) 



where nx is the number density of X and s is the entropy density. Then, when the gaugino 
masses obey the grand-unified-theory (GUT) relation, the yield variable of gravitino is well 
fitted by 



Y 3/2 ~ 2.3 x 1(T 14 x T® 



1 + 0.015 lnTjf } - 0.0009 In 2 T§ ] 
mi/ 2 \ 



.14 v '"i/2 ^(8) 



( 8 ) _l n nnnnir, 2 T( 8 ) 



+ 1.5 x 10~ i4 x — ^ T£> 1 - 0.0371nT^ j + 0.0009 In' T£> , (3.7) 
\m 3 / 2 J L ' J 

where T^ = Tr/10 8 GeV and mi/ 2 is the unified gaugino mass at the GUT scale. The 
above fitting formula agrees well with the numerical result with the error of a few % for 
10 5 GeV <T R < 10 12 GeV. 



#4 See also [53] for possible corrections. 





Case 1 


Case 2 


Case 3 


Case 4 


mi/2 


300 GeV 


600 GeV 


300 GeV 


1200 GeV 


m 


141 GeV 


218 GeV 


2397 GeV 


800 GeV 


A, 














tan/3 


30 


30 


30 


45 


m 


389 GeV 


726 GeV 


231 GeV 


-1315 GeV 


m x\ 


117 GeV 


244 GeV 


116 GeV 


509 GeV 


n (thermal), 2 
1 i LSP (I 


0.111 


0.110 


0.106 


0.111 



Table 1: mSUGRA parameters used in our analysis. 



3.2 Models 



In order to study effects of photo- and hadro-dissociations and p <-» n conversion, it is 
necessary to understand the energy spectra of daughter particles produced by the decay 
of gravitino. If gravitino is not the LSP, gravitino decays into a lighter superparticle and 
standard-model particle(s). Those daughter particles also decay if they are unstable. The 
decay chain strongly depends on the mass spectrum (and coupling constants) of the MSSM 
particles. Since there are so many parameters in the MSSM, it is difficult to consider the 
whole parameter space. Thus, we adopt four parameter points of the so-called mSUGRA 
model to fix the MSSM parameters. Then, for those parameter points, we calculate decay 
chain of the gravitino in details and derive BBN constraints. 

In the mSUGRA model, all the MSSM parameters are determined by the unified gaugino 
mass mi/2, universal scalar mass m , universal coefficient for the tri-linear scalar coupling 
A , ratio of the vacuum expectation values of two Higgs bosons tan/3, and supersymmetric 
Higgs mass /i#. {\hh\ is determined by the condition of electro- weak symmetry breaking.) 
The points we choose are shown in Table [TJ For the points we choose, we also calculate 
the thermal relic density of the lightest neutralino, which is the LSP; the resultant density 
parameters are also shown in Table [TJ (We use h = 0.73 [22J, where h is the expansion 
rate of the universe in units of 100 km/sec/Mpc.) Notice that the Cases 1 and 2 are in the 
so-called "co-annihilation region," the Case 3 is in the "focus-point region," and the Case 4 
is in "Higgs funnel region." 

With the mSUGRA parameters given in Table [JJ, the mass spectrum of the MSSM par- 
ticles as well as coupling parameters are obtained. Then, we calculate the decay rate of 
gravitino r 3 / 2 , taking account of all the relevant decay modes. The lifetime of gravitino is 
given by T3/2 = Tg^; we show the lifetime for each cases in Fig. [TJ The lifetime depends 
on the mass spectrum of the MSSM particles when the gravitino mass is relatively light. 
When gravitino becomes very heavy, on the contrary, r 3 / 2 is determined only by the number 
of possible decay modes. In the present case, as one can see in Fig. [TJ, the lifetime becomes 
insensitive to the mass spectrum of the final-state particles when 777,3/2 ^ 10 TeV. 

The decay processes of the daughter particles from the gravitino decay are also important; 
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Figure 1: Lifetime of gravitino for Cases 1 — 4. The horizontal axis is the gravitino mass. 



we also calculate partial decay rates of the unstable superparticles to follow the decay chain. 
The energy spectra of the decay products are obtained by using Monte-Carlo analysis. At 
the parton level, the decay chain induced by the gravitino decay is systematically followed 
by using ISAJET/ISASUSY packages [23] , while the hadronization processes of partons are 
studied by using PYTHIA package [21]. Then, we obtain energy spectra of (quasi) stable 
particles for the study of the light-element abundances; once the spectra of hadronic particles 
(in particular, those of proton, neutron, and pions) as well as the averaged visible energy 
from the gravitino decay are obtained, we calculate light-element abundances taking account 
of photo- and hadro-dissociation processes as well as p *-* n conversion processes. 



3.3 Numerical results 

We show our numerical results in Figs. [2] — [5]for Cases 1 — 4, respectively. These figures show 
upper bounds on the reheating temperature as functions of the gravitino mass. Since, in this 
section, we consider unstable gravitino, we shade the region where the gravitino becomes 
lighter than the MSSM-LSP and hence is stable. Solid lines show upper bounds on the 
reheating temperature from individual light elements; the upper-left regions are excluded at 
95 % C.L. Compared with the previous work [7], we have newly included contributions of the 
production of the longitudinal component of gravitino. (See Eq. ( 13. 7p ). This modification 
gives us severer upper bounds on the reheating temperature when the gravitino is relatively 
light. Upper bound on the reheating temperature is summarized in Table [2] for several values 
of the gravitino mass . 

When the gravitino is unstable, the LSP (i.e., the lightest neutralino in this case) is 
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Figure 2: BBN constraints for the Case 1 at 95 % CL. Each solid line shows upper bound 
on the reheating temperature from D, 3 He, 4 He, 6 Li, or 7 Li. The dotted line is the upper 
bound on the reheating temperature from the overclosure of the universe. 
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Figure 3: BBN constraints for the Case 2. 
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Figure 4: BBN constraints for the Case 3. 
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Figure 5: BBN constraints for the Case 4. 
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fn 3/2 


Case 1 


Case 2 


Case 3 


Case 4 


300 GeV 


1 x 10 6 ( 3 He) 


4 x 10 5 ( 3 He) 


1 x 10 6 ( 3 He) 


— 


1 TeV 


5 x 10 5 ( 6 Li) 


9 x 10 5 ( 6 Li) 


3 x 10 5 ( 6 Li) 


3 x 10 6 ( 6 Li) 


3TeV 


5 x 10 5 (D) 


4 x 10 5 (D) 


2 x 10 5 (D) 


5 x 10 5 (D) 


10 TeV 


2 x 10 9 ( 4 He) 


2 x 10 9 ( 4 He) 


2 x 10 9 ( 4 He) 


2 x 10 9 ( 4 He) 


30 TeV 


9 x 10 9 ( 4 He) 


8 x 10 9 ( 4 He) 


7 x 10 9 ( 4 He) 


8 x 10 9 ( 4 He) 



Table 2: Upper bound on the reheating temperature (in units of GeV) from BBN for Cases 1 
— 4. The light element which gives the most stringent bound is indicated in the parenthesis. 



produced by the decay of the gravitino. Since the decay of the gravitino occurs after the 
freeze-out epoch of the lightest neutralino, the LSP produced by the gravitino decay survives 
until today. Thus, in this case, we also obtain the upper bound on the reheating temperature 
from the overclosure of the universe. The lightest neutralino has two origins; one is thermal 
relic and the other is non-thermal one from the gravitino decay. However, since the density 
of the thermal relic strongly depends on the MSSM parameters, we do not take into account 
its effect in the calculation of the density parameter. Then, the density parameter of the 
LSP is proportional to Y3/2 and is given by 



V 100 GeV 



n LSP h 2 ~ 2.8 x 10 10 x Y 3/2 ( inn *^ ) , (3i 



where m x o is the mass of the lightest neutralino (i.e., the LSP) which is given in Tabled] 
for each cases. We require that ^lsp^ 2 be smaller than the observed dark matter density: 
^lsp^ 2 < 0.118 (95 % C.L.) [22], and derive upper bound on the reheating temperature. 
The bound is shown in the dotted line in the figures. As one can see, the constraint from 
the relic density of the LSP is less severe than those from BBN unless the gravitino mass is 
extremely large. 

In summary, when the gravitino mass is close the mass of the LSP, the most stringent 
bound is from the overproduction of 3 He, resulting in the upper bound of O(10 5-6 GeV). 
When m.3/2 is around a few TeV, the hadro-dissociation processes become most effective and 
the severest bound is from D or 6 Li. With higher gravitino mass, the constraint is drastically 
relaxed since the gravitino lifetime becomes shorter. In particular, when the gravitino mass 
is larger than ~ 30 TeV or so, the reheating temperature can be as high as ~ 10 10 GeV. 

4 Stable Gravitino 

4.1 General remarks 

Next, we consider the case where gravitino is the LSP and hence is stable. In this case, the 
MSSM-LSP becomes unstable. We assume that the MSSM-LSP is the NLSP and that it 
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decays only into gravitino and standard-model particle(s). Then, the lifetime of the MSSM- 
LSP may become longer than 1 sec. If so, the decay products of the MSSM-LSP induce 
dissociation processes and change the light-element abundances. Since the lifetime of the 
NLSP in this case is approximately proportional to rn\, 2 , the scenario is constrained from 
BBN unless the gravitino mass is small enough. Importantly, the lifetime of the NLSP as 
well as the spectra of charged particles (including hadrons) produced by the NLSP decay 
are strongly dependent on what the NLSP is, so are the BBN constraints. Thus, in the 
following, we discuss several possible cases separately. 

4.2 Case with neutralino NLSP 

First, we consider the case where the NLSP is the lightest neutralino Xv Even in this case, 
in fact, properties of the NLSP depends on MSSM parameters since the lightest neutralino 
is a linear combination of Bino B, neutral Wino W°, and Higgsinos: 

x o = Us b + U W W° + U Su H° u + U & Hl (4.1) 

where H® and H% are up- and down-type Higgsinos, respectively. To make our discussion 
simple, we mostly consider the case that the lightest neutralino is Bino-like (i.e., Ug ~ 1 
and Uyy ~ Ujj u ~ Ujj ~ 0). Notice that the scenarios with Bino-like lightest neutralino 
is realized easily. In particular, if gaugino masses obey the GUT relation, the Bino mass 
becomes lighter than Wino and gluino masses. In addition, assuming radiative electro-weak 
symmetry breaking, the Higgsino mass jjlh becomes larger than M\ in large fraction of the 
parameter space [M]- With those conditions, the lightest neutralino becomes Bino-like. 

If Xi — B, the relevant decay processes of the NLSP is B — > ip^j and B —>■ ip^ff, where 
ipp and / denote gravitino and standard- model fermions, respectively. (See Fig. [6j) Notice 
that the second process is mediated by virtual and on-shell Z-boson as well as virtual photon 
and hence, in our procedure, the decay process B — > ip^Z is taken into account in the process 
B — ► ip^ff. We calculate the decay rate of these processes as well as the energy spectra of 
the final-state particles, in particular, those of partons and electro-magnetic particles. Total 
decay rate of Bino-like neutralino, whose mass is rug ~ Mi, is well approximated by 

r -i ^ T (B -» VV7) + F(S -> Wn (4-2) 

where [551 



T ^^1Y<) = ^^(l-4,) 3 (l + 34 2 ), (4.3) 



'3/2 



sin 0w/3b-^„z ™ 5 

l 3/2 



t(b^^z) = — , _;,r^ -^(i-^ /2 ) 2 (i+3^ /2 ) 



48ttM2 m 



"3/2) K L T °^3/2J 



-X, 



{3 + x| /2 (-12 + x 3/2 ) + x% - x|(3 - x\ l2 )}\ , (4.4) 
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Figure 6: Feynman diagrams relevant for the decay processes of Bino-like neutralino. 



with #w being the Weinberg angle, x 3 / 2 
nig > m 3/2 + m z , 



m 3 / 2 /mg, and Xz = mz/m^. In addition, for 



P 



B^))u.Z 



otherwise. 



l-2fx 



3/2 



+ 4) 



+ {x 



3/2 



X 



2 \2 



1/2 



(4.5) 



while /3s_^ z 

As well as the decay rate and the energy spectra of decay products, it is also necessary to 
obtain the primordial abundance of the lightest neutralino in deriving the BBN constraints. 
The primordial abundance is determined at the time of its freeze-out if there is no entropy 
production after freeze out, and is strongly dependent on MSSM parameters. In our study, 
instead of performing a precise calculation of the primordial abundance assuming a detailed 
model, we first adopt approximated formulas to derive constraints. Following [55J, for the 
abundance in the bulk region, we use 



Y, 



B 



4 x KT 12 x ' m z 



bulk, 



VlOOGeV, 
while, for that in the focus-point and co-annihilation regions, we use 



y B 



9x 10 



-13 



X 



til 



B 



100 GeV 



: focus / co-annihilation. 



(4.6) 



(4.7) 



Constraints on the m 3 / 2 vs. m§ plane with the relic abundance given in Eqs. (J4.6P and 



(14. 7p are shown in Figs. [7J and El respectively. The most stringent constraint is from the 
overproduction of D and 4 He; overproduction of D is due to the hadro-dissociation of 4 He 
produced by the SBBN reactions, while that of 4 He is due to the p <-» n conversion. Compar- 
ing Figs. [7] and [HJ we can see that the constraints from the overproduction of D on the m 3 /2 
vs. nig plane are not so sensitive to Yg. This is because the upper bound on the primordial 
abundance of B from the overproduction of D significantly changes with a sight change of 
the lifetime. (See Figs. O and (TDJ) On the contrary, constraint from 4 He is sensitive to Y§; in 
the focus-point/co-annihilation case where Yg is relatively small, constraint from 4 He is not 
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Figure 7: BBN constraints for the Bino-NLSP scenario with the abundance given in Eq. 
(14. 6p . In the shaded region, gravitino becomes heavier than the lightest neutralino and 
hence is not the LSP. The region with small gravitino mass is allowed. 
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Figure 8: Same as Fig. [3, except for the abundance given in Eq. ( 14. 7ft . 
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important while, in the bulk case, 4 He is too much produced when the Bino mass becomes 
large. 

Next, we treat Y B as a free parameter and derive BBN constraints on 7713/2 vs. rrigYg 
plane. The results are shown in Figs. [9] and \IU\ for m B = 100 GeV and m B = 300 GeV, 
respectively. When the gravitino mass is close to m B , the lifetime of the lightest neutralino 
becomes relatively long. In such a case, the most stringent bound on Y B is from the overpro- 
duction of 3 He due to the photo-dissociation of 4 He. On the contrary, with smaller gravitino 
mass, hadro-dissociation or p <-> n conversion becomes more important than the photo- 
dissociation. In addition, if the gravitino mass is small enough, t b becomes shorter than 
~ 1 sec and, in such a case, decay of the lightest neutralino becomes harmless to the BBN 
scenario. 

So far, we have considered the case where the lightest neutralino is well approximated 
by the Bino. We also considered the case where the lightest neutralino has sizable Higgsino 

0.1, and found that the 



components. We derived constraints for the case Ujj 
constraints are almost unchanged compared to the pure-Bino case. 



U H d ^ 



4.3 Case with stau NLSP 

Next, we consider the case where the NLSP is charged slepton. The basic procedure to 
derive constraints is explained in [TJ]. Here, we concentrate on the case where the NLSP is 
stau f , since, in some models of supersymmetry breaking, like gauge-mediation model [56] 
and mSUGRA model, one of the stau becomes lighter than other sleptons because of the 



renormalization group effects r 5 



In considering long-lived stau it is crucial to take account of the bound state effect which 
significantly change BBN reaction rates, in particular, production rate of 6 Li [20, EHJ [53, dU 
1581 1591 PJ3] . In the previous study by three of the present authors [H] , the number density of 
( 4 Hef ~) bound state was obtained with the use of Saha equation |60j. After the study, there 
have been several works which calculate the number density of bound state by solving the 
Boltzmann equation, which gives more accurate estimate [42j [59[ [43] . In our current study, 
we also solve the full Boltzmann equation to evaluate the number density of bound-state 
taking effect of bound-state (Hf~) into account F 6 



The Boltzmann equations that govern the evolution of the number densities 774 and n\ 
of the bound states ( 4 Hef~) and (Hf~) are 



~dt 

dn\ 
~dT 



-3Hri4 — Tfn 4 + {a r ^v) 
-3Hrii — Tf-rii + {cr r \v) 



\ruu e 



rujrif 



[rip -ntjUf 



/m4 He m f T\ 3/2 _ E , JT 
m p m f T \ e _ Ebl/Tn ^ 

271771! / 



(4.8) 
(4.9) 



# 5 We have checked that the constraints are almost unchanged when the NLSP is selectron or smuon. 
# 6 See also various attempts to tackle this 6 Li problem [52 [6TJ |40l [41] [62J, [44] and a possibility of non- 
standard 9 Be production [63] in the catalyzed BBN induced by the stau. 
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Figure 9: BBN constraints for the Bino-NLSP scenario on the m 3 / 2 vs. m^Y^ plane. The 
Bino mass is taken to be 100 GeV. 
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Figure 11: Time evolution of fractions of r which form bound-states ( 4 Her ) (thick solid 
curve) and (Hf~) (thin solid curve). We also show n 4 /n4 He (thick dashed curve) and n\/n p 



(thin dashed curve). We take Y? = 2 X 10 



-12 



m~ 



1 TeV and I> = 0. 



where n f - is the number density of the free negative charged stau, IV is the decay rate of 
stau, Em — 337.33 keV and Ebi — 24.97 keV are the binding energies of ( 4 Hef ~) and (Hf _ ) 
[58] , and n4 He and n p are the number densities of 4 He and proton including both free and 
bound states. Furthermore, m 4 = m4 He + rrif — Ey± and m\ = m p + rrif — E^i, and the 
thermally averaged recombination cross sections for 4 He and p bound-states are given by 



(a rA v) 
(a r iv) 



98.46 



aE, 



61 



24.62 



m? He v'm4 He T' 
aE bl 



m 2 p Jm p T 



(4.10) 
(4.11) 



An example of the evolution of the abundances of the bound-states ( 4 Hef~) and (Hf~) are 
shown in Fig. [TTJ 

First, we use the following thermal relic abundance of f 



Y f 



7 x 10~ 14 x 



rrif 



100 GeV 



(4.12) 



and derive constraints on m 3/ / 2 vs. vrif plane. The result is shown in Fig. [T2J As one can see, 
the most stringent bound is from the overproduction of 6 Li, which is due to the catalyzed 
process. On the contrary, the bound from D is less severe compared to the Bino-NLSP case. 
This is because, when the stau is the NLSP, the hadronic branching ratio is several orders 
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Figure 12: BBN constraints for the stau-NLSP scenario with the thermal abundance given 
in Eq. (TCT2]) . 



of magnitude smaller than that in the Bino NLSP case, and hence the hadro-dissociation 
processes are suppressed. 

We also treat the primordial abundance as a free parameter and derive upper bound on 
Yf. The results for rrif = 100 GeV and 300 GeV are shown in Figs. [13] and [TU respectively]^!] 
The 4 He-constraint for smaller gravitino mass comes from the p *-> n conversion induced by 
the charged pions which are mainly produced by the decay of tau lepton [TJ]. On the other 
hand, the D-constraint comes from the destruction due to energetic baryons produced by 
the four-body decay. Thus, for the stau case a simple scaling of the constraints on Bino by 
using the hadronic branching ratio does not work. 

As we mentioned, in the present study, we have solved the full Boltzmann equation to 
calculate the number density of the ( 4 Hef~) and (Hf~) bound states. Then, as discussed 
in [591 1^3] . the number density of ( 4 Hef~) is reduced compared to the result with Saha 
equation, resulting in weaker constraint than the previous study. However, the upper bound 
on Yf from the overproduction of 6 Li is increased by the factor 3 or so. Thus, we conclude 
that the previous study with Saha equation provided a reasonable estimation of the upper 
bound on Yf. Notice that, concerning the constraints on the 7713/2 vs. m? plane, the change 



# 7 The bound states with singly-charged nuclei such as (Hf~), (Df~) or (Tf~) may shield coulomb field 
of the nuclei completely and significantly enhance the further reaction rates for these nuclei. These non- 
standard processes might totally have reduced the 6 Li, 7 Li and 7 Be abundances [3911421143] , However, the 
Bohr radius of those bound states are larger than the typical size of the nuclei (or the square-root of the 
cross sections). Then the multi-body problem might be important for reactions and could not have been 
understood well. Therefore, we did not include those effects in the current study. 
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is negligible]^ 8 

4.4 Case with sneutrino NLSP 

Finally we study the case with sneutrino LSP. Although such a scenario is not conventional, 
it is possible to realize sneutrino MSSM-LSP in some gravity mediation scenario. In this 
paper we update the constraint obtained in (13] using the most recent observational data. 
The yield variable for thermally produced sneutrino is estimated as [61] 



y - 2xl °- 14x (ra)- < 413 > 

where mo is the sneutrino mass. With this yield variable, we derive the constraint on 7713/2 
vs. rraj> plane and the result is shown in Fig. [151 The resultant constraint is almost same 
as that in [13] since the most stringent constraints come from overproduction of D and 6 Li 
whose observational abundances have not changed very much. As expected, the constraints 
are much weaker than the Bino- or stau-NLSP cases because the sneutrino dominantly 
decays into particles whose interactions are weak; v — > ip^v. The BBN constraints are from 
four-body decay modes v — > ip^uff, which have very small branching ratios. 

We also treat the primordial abundance as a free parameter and derive upper bound on 
Yp. The result for m^ = 300 GeV is shown in Fig. [161 Compared to the previous study 
[T3] , it is seen that the constraint from 4 He is slightly milder due to the larger observational 
constraint on Y p than that adopted in |13j . 

5 Conclusions and Discussion 

In this paper we have obtained the improved BBN constraints on both unstable and stable 
gravitino cases. We have taken into account recent theoretical and observational progresses 
in the study of BBN constraints on long-lived unstable particles. 

In the unstable gravitino case, we have included effects of the longitudinal component 
which was not considered in the previous studies, which have led to a more stringent upper- 
bound on the reheating temperature for gravitino with relatively small mass. The upper 
bound on the reheating temperature strongly depends on the gravitino mass and has been 
summarized in Table [2] for the MSSM parameters used in our analysis. We can see that the 
upper bound has mild dependence on the mass spectrum of the MSSM particles, and that the 
bound becomes almost independent of the MSSM parameters when the gravitino becomes 
much heavier than the MSSM superparticles. Our results give very stringent constraint on 
the scenario of thermal leptogenesis [65] which requires the reheating temperature higher 
thanO(10 9 GeV) jgSlEZ]. 



# 8 Ref. [40] suggested a mechanism to reduce the 7 Li ( 7 Be) abundance through ( 7 Bef~) +p — ► ( 8 Bf _ ) +7 
witch is effective for Y? ^ 10 -11 and Vf ^ 10 3 sec. Because we have adopted the conservative observational 
7 Li abundance which agrees with the SBBN prediction, this mechanism does not change our results. 
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given in Eq. (j4.13p . 
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When the gravitino is the LSP and stable, the constraints depend on which superparticle 
is the NLSP. We have considered Bino, stau and sneutrino as NLSP and obtained constraints 
on their properties. Among these possibilities, the hadronic branching ratio becomes largest 
when the Bino is the NLSP, for which the constraints from the hadro-dissociation processes 
are the most stringent. For the case of stau as NLSP, we have estimated the effects of bound- 
state of stau with 4 He accurately by solving the Boltzmann equation describing evolution 
of the bound-state abundance. However, the bound from the catalyzed production process 
of 6 Li is almost unchanged with the previous study where Saha equation was used, and the 
overproduction of 6 Li gives stringent bound on the primordial abundance of stau. In these 
two cases, the constraints are so stringent that the gravitino mass larger than ~ 10 GeV is 
excluded when the mass of the MSSM-LSP (i.e., Bino or stau) is lighter than 1 TeV. If the 
sneutrino is the NLSP, on the contrary, BBN constraints become drastically weaker since the 
sneutrino mainly decays into very weakly interacting particles, i.e., gravitino and neutrino. 
However, in any case, our results provide severe constraint on some scenario of gravitino 
dark matter. Indeed, as we can see in Figs. [7J [HI d21 and (TBI the constraints from BBN are 
severer than that from the overclosure of the universe. This fact implies that the abundance 
of the gravitino produced by the decay of the MSSM-LSP cannot be sufficient to realize the 
gravitino dark matter if the thermal relic abundance of the MSSM-LSP is assumed. Thus, 
in order to realize the gravitino dark matter, most of the gravitino should be produced by, 
for example, the scattering processes of thermal particles at the time of the reheating after 
inflation [9] or by the decay of scalar condensation [331 150] . 
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